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ABSTRACT 

Aims. In the frame of the search for extrasolar planets and brown dwarfs around early-type main-sequence stars, we present the results 
obtained on the early F-type star 6>Cygni. 

Methods. ELODIE and SOPHIE at Observatoire de Haute-Provence (OHP) were used to obtain 91 and 162 spectra, respectively. 
Our dedicated radial-velocity measurement method was used to monitor the star's radial velocities over five years. We also use 
complementary, high angular resolution and high-contrast images taken with PUEO at CFHT. 

Results. We show that #Cygni radial velocities are quasi-periodically variable, with a ^ 150-day period. These variations are not due 
to the ^ 0.35-Mq stellar companion that we detected in imaging at more than 46 AU from the star. 

The absence of correlation between the bisector velocity span variations and the radial velocity variations for this 7 km s v sin i star, 
as well as other criteria indicate that the observed radial velocity variations are not due to stellar spots. The observed amplitude of 
the bisector velocity span variations also seems to rule out stellar pulsations. However, we observe a peak in the bisector velocity 
span periodogram at the same period as the one found in the radial velocity periodogram, which indicates a probable link between 
these radial velocity variations and the low amplitude lineshape variations which are of stellar origin. Long-period variations are not 
expected from this type of star to our knowledge. If a stellar origin (hence of new type) was to be confirmed for these long-period 
radial velocity variations, this would have several consequences on the search for planets around main-sequence stars, both in terms 
of observational strategy and data analysis. 

An alternative explanation for these variable radial velocities is the presence of at least one planet of a few Jupiter masses orbiting at 
less than 1 AU; however this planet alone does not explain all observed features, and the (fCygni system is obviously more complex 
than a planetary system with 1 or 2 planets. 

Conclusions. The available data do not allow to distinguish between these two possible origins. A vigourous follow-up in spectroscopy 
and photometry is needed to get a comprehensive view of the star intrinsic variability and/or its surrounding planetary system. 

Key words, techniques: radial velocities - stars: early-type - stars: planetary systems - stars: individual: #Cygni 



* Based on observations made with the ELODIE and SOPHIE spec- 
trographs at the Observatoire de Haute-Provence (CNRS, France) and 
with the PUEO adaptive optics system at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
of Canada, the Institut National des Sciences de l'Univers of the Centre 
National de la Recherche Scientifique of France, and the University of 
Hawaii. 

** Tables of radial velocities are only available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strabg.fr/cgi-bin/qcat ?J/A+ A/ 



1. Introduction 

Radial-velocity (RV) surveys have lead to the detection of more 
than 300 planets during the past decad^J These surveys mainly 
focus on solar and later-type main-sequence (hereafter MS) 
stars (^ F7) which exhibit numerous lines with low rotational 
broadening, making them ideal targets for classical velocimetry. 
However, it is crucial to understand how planetary systems form 
over a wide variety of parent stars, and to know in particular 
if there is a correlation between the planet masses and the par- 
ent star masses as predicted for instance by Kennedy & K enyon| 
(2008), and to constrain formation models such as those from 



1 A comprehensive list of known exoplanets is available at 
http://exoplanet.eu 



2 



M. Desort et al.: Extrasolar planets and brown dwarfs around A-F type stars 



Ida & Lin (2005) (but also Hop 2006 and |LaugMin"etaLl 2004 
for M dwarfs), and/or a correlation between the planet occurence 
and the parent star masses. Concerning massive stars, surveys of 
subgiant/giant stars have started to provide first information on 
planets at orbital distances typically greater than 0.7 AU (e.g., 
IJohnson et all 20 06; IJohnson et all 200 7; IHatzes et al.l 2005; 
INiedzielski e HD2007; [Lovis & Mayor| 2007; lSato_e£jD2008). 
Closer separations have to be investigated by observing massive 
main-sequence stars. In this frame, we developed a tool dedi- 
cated to the search for planets around early (A-F) type stars. The 
method allowing the measurement of the RV of rapid rotators is 
described by Galla nd et al.l (2005a; hereafter Paper I). 

We started in 2005 two surveys dedicated to the search for 
extrasolar planets and brown dwarfs around a volume-limited 
sample of A-F main-sequence stars i) with the ELODIE fiber- 
fed echelle spectrograph (Bara nne et al.l 1996) mounted on 
the 1.93-m telescope at the Observatoire de Haute-Provence 
(OHP, France) in the northern hemisphere, and ii) with the 
HARPS spectrograph ( |Pepe et aT] 2002) installed on the 3.6-m 
ESO telescope at La Silla Observatory (Chile) in the southern 
hemisphere. In 2006, the ELODIE spectrograph was replaced by 
SOPHIE ( |Bouchy et aLl2006). We d etected with ELODIE a planet 
around an F6V star (IGalland et al.l 2005b, Paper II) and a brown 
dwarf around an A9V star (IGalland et al.l 2006b, Paper IV), 
and with HARPS a two-planet system around an F6IV-V star 
dDesort et al.1 2008. Paper V). We also derived the first statistics 
of planet existence around A-F stars thanks to our HARPS sur- 
vey ( |Lagrange et al.| 2009, Paper VI). 

We present and analyse in this paper the RV variations of 
8 Cygni. Section|2]provides the stellar properties, and the various 
data obtained on this object. In Sect. [3] we discuss the origin of 
the observed RV variations. 



2. Stellar characteristics and measurements 

2.1. Stellar properties 

flCygni (HD 185395, HIP 96441, HR7469) is a Mi = 
1.38 + 0.05 M star, with an age estimated to 1.5^ Gyr 
dNordstrom et al.12004), and lo cated at 18.33 ± 0.05 pc from the 
Sun dESAI 1997. Ivan Leeuwenl 2007). We took its rotational ve- 
locity v sin i, effective temperature T e g, and surface gravity log g 
from |Erspamer & North (2003) and Gray et al. (2003) (values in 
Table [1J). We assume a spectral type F4V, commonly attributed 
to this star as, e.g., in the Bright Star Catalogue (Hoffiei t~et al.l 
1991) or in the HIPPARCOS catalogue (ESA]1997). 



Table 1. 8 Cygni stellar properties. Photometric and astro- 
metric data are extracted from the HIPPARCOS catalogue 
(ESAj 1997, Ivan Leeuwenl 2007); spectroscopic data are from 
INordstrom et~aT1 (2004) and |Erspamer & North| (2003). 



2.2. Spectroscopic data 
2.2.1 . Description of the data 

Between 2003 and 2006, we recorded 91 high S/N spectra 
of 8 Cygni with ELODIE and, between November 2006 and 
December 2008, we recorded 162 spectra with SOPHIE. 

The wavelength range is 3850-6800 A for ELODIE and 
3872-6943 A for SOPHIE. Typical exposure times were 15 and 
3 min respectively for ELODIE and SOPHIE, leading to a signal- 
to-noise ratio S/N ~200. The exposures were performed with 
simultaneous-thorium spectra to follow and correct for the possi- 
ble drift of the instrument due to local temperature/pressure vari- 
ations (whose impact shows a standard deviation of 2.5 ms -1 ). 
With SOPHIE we used the high-resolution (R * 75 000) mode. 



Parameter 




6 Cygni 


Spectral Type 




F4V 


v sin i 


[kms- 1 ] 


7 


V 




4.49 


B-V 




0.395 ± 0.015 


n 


[mas] 


54.54 + 0.15 


Distance 


[pc] 


18.33 ±0.05 


M v 




3.14 


[Fe/H] 




-0.08 


Teff 


[K] 


6745 


l°gg 




4.2 


Mi 


[M Q ] 


1.38 ± 0.05 


Age 


[Gyr] 


1 5 +0 - 6 


ppm (a) 1 ' 


[mas yr -1 ] 


-8.15 


ppm ((5)+ 


[mas yr -1 ] 


-262.99 



' the proper motion is affected from the orbital motion that we discuss 
in Sectionl2~3l 



2.2.2. Radial velocity variations 

The radial velocities (Fig.[T]i are measured using a dedicated tool 
(S AFIR) described in Paper I and based on the Fourier interspec- 
trum method developed in lChellil (2000). The uncertainty asso- 
ciated with ELODIE data is 9ms~ 1 on average, consistent with 
the value obtained from our simulations (see Paper I). In the case 
of SOPHIE data, the uncertainty is 5 m s on average (taking the 
photon noise and instrument stability into account). 

As 8 Cygni has a relatively low projected rotational-velocity 
(v sin i = 7 km s ), we could also measure the RV using a gaus- 
sian adjustement to the cross-correlation function (CCF). The 
results obtained by the two different methods are found to be 
consistent. 

ELODIE and SOPHIE RV data show (Fig.[T]i a quasi-periodic 
signal with peak-to-peak amplitude of about 220 ms -1 , much 
larger than the uncertainties. A drift in the RV curve seems more- 
over to be present over the whole data set, which could be at- 
tributed to a stellar companion (Section |2~3j ). We finally note that 
the amplitude of the RV variations could in addition be slightly 
variable. 

We used the CLEAN algorithm (Hog bom 1974), applied to 
Lomb-Scargle periodograms to derive the periodograms of the 
radial velocities measured with both instruments (Fig. [2}; this 
algorithm removes the aliases associated with temporal sampling 
of the data: it deconvolves the window function iteratively from 
the initial "dirty" periodogram to produce the resulting cleaned 
periodogram. 

In the case of ELODIE data, the periodogram shows one peak 
at a period of 128 ± 5 days (the uncertainty is evaluated with the 
full width at half maximum of the highest peak). In the case of 
SOPHIE data, the peak corresponds to a period of 158 + 10 days. 



2.2.3. Line profile variations 

SOPHIE Lomb-Scargle periodograms of the BVSs and curva- 
tures (defined as in IHatzesI 1996) are presented in Fig. [3] to- 
gether with false alarm probabilities (FAP, Kiir ster et al.l 1997). 
In the case of the SOPHIE data, a peak is seen at approximately 
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Fig. 1. Radial velocities of 6»Cygni obtained with ELODIE (left) and SOPHIE (right). 
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Fig. 2. ELODIE (top) and SOPHIE (bottom) CLEAN peri- 
odograms of the radial velocities. On each data set one peak is 
observed, either at 128 + 5 days or at 158 ± 10 days respectively. 



140 days, i.e., not very different from the one measured in the 
SOPHIE RV variations. No peak is detected on the yet noisier 
ELODIE data. 

Figure |4] presents the bisector velocity spans (BVSs) as a 
function of RVs for ELODIE and SOPHIE spectra. The amplitude 
of the B VS variations is quite small: 50 m s _I , much smaller than 
the amplitude of the RV variations. No clear correlation is seem 
between the BVS and RV variations. It shows that the spectra 
are mainly shifted in radial velocity without significant changes 
in the lines shape. 
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Fig. 3. ELODIE and SOPHIE Lomb-Scargle periodograms 
of the bisector velocity spans (top 2 panels) and curva- 
tures (bottom 2 panels) with false alarm probabilities. On 
SOPHIEperiodograms the period of 158±10 days is repre- 
sented with vertical dashed lines (dotted for the ±10 range). 



2.2.4. Stellar jitter 

Finally, to quantify its short-term variations, we monitored 
#Cygni for 1 hour in a row (as we did for HD 60532, Paper V). 
High-frequency variations are due to stellar phenomena and pro- 
duce a noise (jitter) that has to be taken into account in the anal- 
ysis of longer period variations. In October 2008, 46 consecu- 



tive spectra were then taken with SOPHIE under average obser- 
vation conditions (airmass below 1.2, S/N - 160). The result- 
ing RVs and associated bisector velocity spans are presented in 
Figure[5] It appears that the short-term variations can account for 
an RV amplitude of ^30 ms -1 (cr rv = 6.4 ms -1 ), and that the to- 
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Fig. 4. Bisector velocity span versus RVs for ELODIE (top) and 
SOPHIE (bottom) data, showing that there is no correlation be- 
tween line profiles and radial velocities. The dashed lines on 
SOPHIE data show the minimum effect of short-term variations, 
see Section l2.2.4l 



tal bisector velocity span amplitude over the whole SOPHIE data 
(Fig. [4] bottom) set can be explained only with those short-term 
variations. In the following attempt to analyse the high ampli- 
tude RV variations, we will adopt increased RV uncertainties (at 
least ±6.4 ms , fixed to that value) in order to take this stellar 
jitter into account. 



2.3. Imaging data 

We observed flCygni at high angular resolution and high 
contrast with the adaptive-optics (AO) instrument PUEO 
(Rigautetal. 1998) mounted on the 3.6-m Canada-France- 
Hawaii Telescope (CFHT, USA). The near-infrared camera used, 
KIR ( |Doyon eFaT1 l998), has a field of view of ~ 35" x35", with 
a seal 



e of ~ 0.035" per pixel. We performed non-saturated ex- 
posures as well as 30-second saturated ones, to investigate at a 
deeper level the vicinity of the star. A special care was taken to 
ensure that the non-saturated exposures could be used as refer- 
ences for an accurate measurement of the possible companions 
positions relative to the central star and as references for the 
measurement of the photometric contrast between the star and 
the possible companions. 

SCygni was observed in June 2004, September 2005 and 
November 2007. The log of observations is given in Table [2] 
A classical reduction has been performed using the software 
ECLIPSE dDevillardl l997). A candidate companion (CC) is seen 
in non-saturated images (see Figure [6} with the narrow band- 
width Fen filter (A - 1.644//m, AA - 0.015 yum). A deconvolu- 
tion algorithm using the method described in lVeranl (1997) has 
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Fig. 5. Top: RVs for 1-hour SOPHIE data. Bottom: Bisector ve- 
locity spans versus RVs for 1-hour SOPHIE data. 

Table 2. Log of WDS observations between 1889 and 1968 fol- 
lowed by the AO observations using PUEO at CFHT between 
June 2004 and November 2007. In the case of WDS data the 
uncertainties are unknown. 



Date 


P 


e 


Contrast 


Band 




["] 


[°] 


[mag] 




1889.37 


3.62 


43.9 






1892.38 


3.79 


47.0 






1898.46 


3.37 


49.2 






1898.63 


3.71 


46.9 






1958.58 


3.42 


51.7 






1968.72 


2.92 


59.9 






2004-06-28 


2.510 ±0.021 


67.37 ± 0.48 


4.6 ±0.1 


Fen 


2007-11-16 


2.369 ± 0.005 


69.02 ±0.11 


4.6 ±0.1 


Fen 



been applied to derive the contrast and angular separation (p) 
between the star and the companion. 

Fig. [7] shows the relative positions of the CC between 2004 
and 2007. Clearly, the CC is not a background star but it is bound 
to <9Cygni. Given the HIPPARCOS distance, 18.6 pc, we derive a 
projected separation of 46.5 AU between the two objects, thus a 
minimum period of roughly 230 years, assuming a circular orbit. 
Moreover, we see in Fig. [7] that the orbit of the companion is 
— still very partially — resolved over a three year period of 
observation. 

The measured contrast between #Cygni and its companion 
#Cygni B is 4.6 + 0.1 mag in H band; in K band, the mea- 
sured contrast is 4.5 + 0.1 mag. Given the star apparent mag- 
nitudes provided by Skruts kie et al.l (2006), and distance (see 
above), we deduce H and K absolute magnitudes of 7.0 + 0.1 and 
6.7 ±0.1 mag respectively. Using the BCA98 evolutionary mod- 
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Fig. 6. High contrast and angular resolution image acquired with 
the adaptive optics system PUEO installed at CFHT: a compan- 
ion is clearly visible. 



1000 - 



800 - 



to 600 - 
< 

400 



200 




O liiiiliiiiliiiiliiiilii.il 
-2150 -2200 -2250 -2300 -2350 -2400 

Aq (mas) 

Fig. 7. Temporal evolution of the separation of the stellar com- 
panion of 8 Cygni with respect to the central star. Filled symbols 
represent the position of the stellar companion relative to the 
main star for the three AO observations (filled circle: June 2004; 
filled square: September 2005; filled triangle: November 2007), 
whereas empty symbols show the positions at the same dates in 
case of a background star. The curve shows the path that would 
have followed the candidate companion if it was a background 
star, taking the star proper- and parallactic -motion into account. 



els ((Baraffe et al. 1998), and assuming any age above ~100Myr, 
we deduce a mass ma — 0.35 M0 for the companion (the evolu- 
tionary effects are negligible). Using the empirical relation given 
by Delfos se"et al.l (2000) for Mh — 7, we find a comparable mass 
ni2 - 0.33 M0 for the companion. 

8 Cygni was classified as a double star ( Dommang et & Nys| 
1994). The Was hington Double Star (WDS) data 
( |Hartk opf & Mason] 2001) indicate a visual companion de- 
tected several times since 1889, with a magnitude of ~12, i.e., 
comparable to the visual magnitude expected from a 0.35- 
M0 star at 8 Cygni distance. The relative position of this object 



1 0-1-2 -3 -4 

Act (arcsec) 

Fig. 8. Positions of the companion star with respect to the pri- 
mary star for more than one century. The star represents the pri- 
mary 8 Cygni, the circles the WDS data (uncertainties are un- 
known), and the triangles our AO data (uncertainties are smaller 
than the symbols). 



varies between 1892 and 1968 (see Tab.[2]i, and indicates that 
this companion is bound to 8 Cygni. The companion was not 
detected by HIPPARCOS because the contrast with 8 Cygni was 
too high. It is reasonable to say that the WDS companion and 
the one found with PUEO is the same. Then we can track its 
motion over more than a century (Fig. [8}, but its orbit is still 
very incomplete. 

3. Origins of the observed RV variations 

We investigate hereafter different possible origins for the RV 
variations: stellar phenomenon (spots, pulsations), and planets. 
Beforehand, we estimate the possible impact of the stellar com- 
panion on the spectroscopic data. 

3. 1 . Impact of the stellar companion on the RV data 

8 Cygni B could a priori play a role on the measured radial veloc- 
ities in two ways. Firstly, given their separation (-2"), and given 
the usual seeings and the entrance width of the optical fiber, 
the spectra of the two stars are superimposed in the ELODIE or 
SOPHIE data (fiber entrance of 2" and 3" respectively). This 
could introduce a bias in the measurements. However, the con- 
trast of 4.6 magnitude in H band translates to a contrast of 
7.9 mag (a flux ratio of ~ 1500) in V band. Hence, the signal 
of the secondary is negligible in our spectra (we expect a radial- 
velocity effect below 1 ms _1 ). If very actuve, it could still pro- 
duce a weak Ha signature superimposed on the 8 Cygni Ha line. 
This does not affect our results as this line is not taken into ac- 
count for the RV measurements. 

With the classical cross-correlation technique, the potential 
pollution of the spectrum by a stellar companion can be tested 
using various masks. We therefore checked that the RV ampli- 
tudes remain identical when using various masks. Also, the RVs 
are identical if we use either the red part or the blue part of the 
spectra to measure them. This confirms that the spectrum of the 
companion has no impact on the measured radial velocities. 



6 



M. Desort et al.: Extrasolar planets and brown dwarfs around A-F type stars 




Unprojected separation r [AU] 

Fig. 9. Possible drift induced by the stellar companion given the 
unprojected true separation r to the primary star. 



Secondly, the stellar companion induces of course radial ve- 
locity variations of the primary. Given the companion properties, 
and assuming the system seen edge-on, we estimate the maxi- 
mum drift possibly induced on the primary star by plotting the 
acceleration projected on the line of sight z, with respect to the 
true separation r (Fig. 19b 



dh 
dt 2 



Gm 



2 r ■ (P\\ 

cos y arcsin \—)\, 



were G is the gravitational constant, ni2 — 0.35 Mq the mass 
of the secondary, and p = 46 AU the closest projected separa- 
tion measured. We find a maximum drift of 12 m s _1 yr _1 , which 
would lead to a maximum drift of 60 ms -1 over 5 years. In fact, 
a similar drift has to be included into our fit of the radial-velocity 
curve (see Section pO) . SCygni B could then explain such a 
drift, but of course the observed periodic radial-velocity varia- 
tions with an amplitude larger than 150 ms -1 are not explained 
by the presence of this stellar companion. 



3.2. Stellar phenomenon 
3.2.1. Stellar spots 

We saw in the previous section that there is no correlation be- 
tween the star RV variations and the BVS variations. Given the 
star projected rotational-velocity, the instrument resolution, and 
according to the study presented in lDesort et al.l (2007). we can 
definitely conclude that the observed RV variations are not due 
to stellar spots. Indeed, would this be the case, a correlation be- 
tween the BVS and the RV would be observed. Typically, given 
the star's properties, one or more spots on a inclined star would 
be needed to reproduce a periodic signal, and a linear correla- 
tion coefficient -0.5 between the bisector velocity spans and 
the RVs. Given the observed RV amplitude, the amplitude of the 
BVS would therefore be much higher than what is actually ob- 
served. 

Also, in such a case, one would expect significant photo- 
metric variations. A single spot producing such an RV variation 
would induce a photometric amplitude between 5 and 30 mmag, 
depending mainly on the star inclination and the spot loca- 
tion dDesort et al.l 2007). The photometry given by HIPPARCOS 
(ESA 1997) is constant with a scatter of only 0.004 mag. We 
recognize however that the HIPPARCOS data were not recorded 
simultaneously with the spectroscopic ones, so this photomet- 
ric argument is certainly weaker than the absence of correlation 
between RV and BVS variations. 

No clear emission in the core of the Ca n lines is observed 




3920 3930 3940 

Wavelength [A] 

Fig. 10. No emission is observed in the Ca n K line for the 
#Cygni spectra. This is the average spectrum of all the 
SOPHIE spectra used, after recentering by the RV variations 
measured. 



We looked for possible long-term stellar variations using 
classical Ha,/?, y indicators. Very faint Ha variations are de- 
tected, but they are not correlated with the RV variations. As 
the Ha indexes are moreover quite sensitive to pollution by the 
thorium lamp, we cannot attribute them a stellar origin. 

Finally, the stellar rotational period is less than 7 days ac- 
cording to its v sin i and assuming a stellar radius typical for this 
type of star. Hence a spot or a set of spots can not explain RV 
variations with periods of one hundred days or more. 

We can therefore safely conclude that spots are most proba- 
bly not responsible for the observed RV variations. 

3.2.2. Granulation 

Another source of variation could be the attenuation or sup- 
pression of the convective blueshift due to the presence of 
plages (e.g., |Deming & Plymate| 1994, |Marquez et al.| 1996). 

The expected convective blueshift expected for 
#Cygni would be in the range 400-800 ms 1 depending 
on the line (see Gray 2009 for a F5IV-V, which is the closest 
to flCygni in his sample). The observed variations could 
then be due to a cyclic variation with plages filling factor 
varying between and 30% typically in order to produce 
the 220 ms 1 variation. However, such a variation of the 
convection properties should lead to a strong variation of 
the bissectors and of the Ca n index, which is not the case. 



(see Fig. 10 for the Can K line); this excludes a high level of 
activity. 



3.2.3. Stellar pulsations 

Generally, pulsations induce line-profile variations which 
strongly affect the bisector velocity span (Paper VI). Besides, 
the timescale of the observed radial-velocity variations (> 100 
days) is far larger than the ones of pulsations known for this type 
of main-sequence stars (it is in fact more characteristic of giant 
stars variability). Variability periods of a few days are observed 
in the case of the pulsating y Doradus stars (Mat hias et al. 1 2004). 
Moreover, if we integrate the RVs measured between a mini- 
mum and a maximum of the amplitude for half a period, we end 
up with a total displacement close to the stellar radius, which, if 
even possible, would lead to detectable photometric variations. 

It is therefore unlikely that classical pulsations are respon- 
sible for the observed RV variations. However, the presence of 
a peak at about 150 days in the BVS periodogram and in the 
RV periodogram indicates that the period of the RV variations is 
linked in some way to the low amplitude line shape variations. 
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Such a situation has never been reported to our knowledge and 
is indeed quite puzzling. We cannot at this stage exclude that we 
could be facing a new type of stellar variability, undetected so 
far because of a lack of long-term, very-precise RV monitoring 
of main-sequence stars. 

3.3. Planet(s) around 6 Cygni 

6 Cygni RV curve shows quasi-periodic variations with a pe- 
riod of ~ 130-150 days, together with a variable positive trend. 
Moreover, the RV curve seems to be modulated in amplitude. 
We now try to see whether these variations can be attributed to a 
planet or a planetary system. 

We first try to fit the whole data set (with uncertainties set 
to +6.4 m s _1 ) assuming a single planet and allowing a drift. We 
fail to reproduce correctly the observed RVs. Figure 1 1 shows an 
example of fit with a planet on a circular orbit (plus a drift). The 
planet that produces such an RV curve has a mass of 2.3 Mj up , 
and is on a 0.6 AU circular orbit (Table|3]l. The algorithm ends up 
with a period of 155 days hung on the SOPHIE data set, but fails 
to fit properly the ELODIE data set, as if there was period/phase 
change with time (clearly visible as structure in the residuals 
near the range [3600-3800] days). Note that the additional drift 
needed is 16 m s yr~' , slightly higher than the maximum value 
that the binary companion would probably produce. We note that 
such a system which fails to fit satisfactorily the whole set of 
data (taken over more than 5 years) would allow to fit the data if 
they were limited to one or two consecutive periods. 

We then tried to fit simultaneously the whole set of data, as- 
suming the presence of several planets, and using a genetic al- 
gorithm search. No stable solution was found. The only satisfac- 
tory fits are achieved by unstable systems with orbits that cross 



HD185395 



ELODIE+SOPHIE 



each other. We show an example in Fig. 12 of a fit obtained with 
a two-companion plus drift model. The residuals are still very 
high and the fit is not improved compared to the one-companion 
plus drift model (residual rms 35 ms -1 versus 39 ms -1 , and the 
drift that we get is approximately the same: 17 ms -1 yr ver- 
sus 16 ms -1 yr -1 ). Moreover, such kinds of configurations with 
massive planets on so close orbits are not dynamically stable. 

The system is then obviously more complex than just con- 
sisting of one, two or even three planets plus a drift. We explored 
then more exotic configurations: 

- Instead of harbouring one single planet, the system could 
consist of a binary planet system orbiting the star, much like 
the Pluto-Charo system, but with larger masses. However, 
the radial-velocity signal generated by this configuration 
would be very close to the one generated by a single planet 
having the total mass. The only departure from the pure 
Keplerian system would be related to the secular perturba- 
tions of the orbit due to the binary nature of the system. The 
associated secular period would be at least several hundreds 
of primary orbital periods, i.e., much longer than our obser- 
vation time span. Hence we should mainly detect the primary 
orbital signal with good accuracy. 

- One could also think of 2 co-orbiting planets locked in 1:1 
mean-motion resonance. Such a configuration has been ob- 
served in the system of Saturn satellites. The two satellites 
orbit Saturn on the same orbit. In the rotating frame, their 
have synchronised horseshoe-like libration motions that pre- 
vent them to collide. The less massive one has the largest am- 
plitude motion ( Yo der et al.l 1983). The associated libration 
period is here again a few hundreds of orbital periods of the 
primaries. So, over a smaller time span, the radial-velocity 
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Fig. 11. ELODIE and SOPHIE radial velocities and orbital solu- 
tion for Cygni, with one planet (2.3Mj up at ~0.6 AU, circular 
orbit) and a drift. The residuals to the fitted orbital solution are 
displayed below. 

signal of the whole system will mainly consist of the sum of 
the individual signals of the two planets. If we assume that 
the two planets have zero eccentricities, the radial velocity 
signal will have the form A cos(nf ) + B cos(nf + <f>), where 
n is the common mean motion of the two planets, A and B 
are amplitudes related to the masses of the planets, and <p is 
a phase shift that depends on their current mutual configura- 
tion. This can be rewritten as C cos(nf +</') where C and ifr are 
new amplitudes and phases that depend on A,B and <p. This is 
equivalent to the signal generated by one single planet. This 
remains true if the planets have small eccentricities; and if 
they have larger eccentricities, the mutual system is not sta- 
ble. Here again, the only departure from this signal will be 
the temporal variation of C and ijj that is related to the mu- 
tual libration motion of the two bodies. Hence we should not 
expect to see changes before several decades. The signal de- 
tected over our observation timespan should not be different 
to that of a single planet. 



4. Concluding remarks 

The radial velocities obtained with ELODIE and SOPHIE on 
6 Cygni are quasi-periodically variable, with a ^ 150 days pe- 
riod. We have investigated different possible origins to these RV 
variations, either star or planet related. 

Given today knowledge on stellar activity, using several cri- 
teria as usually done in this type of studies, we fail to attribute 
these variations to the star itself (spots, pulsations). However the 
presence of a peak in the periodogram of the BVSs at 140 
days, i.e., close to the period of the RV variations prevents us 
from excluding totally a stellar origin. 

We tried to fit the data with one planet orbiting at less than 
1 AU and with a mass of a few Jupiter masses, taking also into 
account the impact of the observed companion star. It appears 
that such an hypothesis allows us to fit only part of the data, 
recorded over a limited number of periods, but fails to fit satis- 
factorily the whole set of data (taken over more than 5 years). 
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Table 3. ELODIE/SOPHIE best orbital solution for Cygni, con- 
sidering one planet and a drift. 



Parameter 




6» Cygni b 


P 


[days] 


154.5 ± 0.4 


T 


[JD-2450000] 


4016 ± 1 


e 




(fixed) 


CO 


[deg] 





K 


r -It 

[ms '] 


70 ± 4 


N 




253 


0"o-c 


[ms" 1 ] 


38.7 


reduced x 2 




6.1 


a\ sin i 


[10~ 3 AU] 


0.99 


f(m) 


[10- 9 Mq] 


5.5 


M, 


[Mq] 


1.38 


mi sin i 


[M Jup ] 


2.29 


a 


[AU] 


0.63 


drift 


[ms -1 yr -1 ] 


16 ±4 
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Fig. 12. ELODIE and SOPHIE radial velocities and orbital so- 
lution for # Cygni, with two planets (2.1 Mj up at ~0.6 AU and 
0.8 Mj up at ~0.5 AU on circular orbits) and a drift. The residuals 
to the fitted orbital solution are displayed below. 



More complex systems were investigated, but no convincing re- 
sult was obtained. More observations and detailed studies of 
gravitational interaction between the 2 planets are needed to un- 
derstand this sytem. 

With the data available, we are then not able to conclude on 
the origin of these puzzling RV variations. 

If a planet origin is confirmed, then 9 Cygni, with a spectral 
type of F4V, would be the earliest main-sequence star hosting 
planets found so far. Moreover, it would be one of the few low 
metallicity stars hosting planets. Its planetary system would not 
be simple, and would in particular include strongly interacting 
planets. 

If a stellar origin was to be confirmed, then this would show 
that, unexpectedly, some main-sequence stars, not classified as 
active from usual criteria {e.g., CanH&K indexes) or pulsating, 
may undergo intrinsic variations which produce quasi periodi- 
cal, large-amplitude and long-period (more than 100 days) RV 
variations, with at the same time, low levels of line shape de- 



formations (hence small amplitude bisector velocity span vari- 
ations). Such situations have not been considered so far in the 
analysis of RV variations, and would need to be considered in 
future searches for such long period planets, both in terms of 
observational strategy and data analysis. 

In any case, we can conclude that 6 Cygni is an individual 
complex system that deserves much more observations to be 
understood, and that it may also serve as an example for other 
searches. 
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